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A vortex theory of screw/propellers with variable cir- 
culation according to the blade and its azimuth is proposed, 
the problem is formulated and circulation is expanded in a 
Fourier series. Equations are given for inductive velocities 
in space for screws, including those with an infinitely large 
number of blades! and expansion of the inductive velocity by 
blade azimuth of a second screw. Multiparameter improper in- 
tegrals are given as a combination of elliptical integrals and 
elementary functions, and it is shown how to reduce elliptical 
integrals of the third kind with a complex parameter to inte- 
grals with a real parameter. 
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List of Symbols/Notations/Abbreviations 


x,y,z — coordinate system 
V,V x ,Vy — flight speed 

R,r — radii of the screw and blade cross section 
ip — blade azimuth 
k — number of blades 
a Q — angle of "blade cone" 

— inductive velocity averaged over the disk of the screw 
r(r,ip), r (r), r nc (r), r ns ( r ) — circulation of airspeed around 

blade cross section and its har- 
monics 

v, v Q (r), v mc ( r )> v ms (r) — inductive velocity of screw F(r,^) and 

its azimuth 0 harmonics 

ft — angular velocity of rotation of screw around its axis 
w — inductive velocity of abstract screw of radius p with constant 
circulation y(ip) over the radius 

W n ,W n’0^ r ^’ w n,mc^ r ^’ W n ms^ — inductive velocity of abstract 

screw with circulation y =e^ n ^ and 

n 

its harmonics 

X 2 ?y 2 ,z 2 — coordinates of center of hub of second screw 
aQ 2 — value of a^ for second screw 

sgnA — unit with sign of value of A, sgnO=l in this case 


F q = dF I dq 

K = dK/dq 

V = ^ ^ = = y - a„p - 

'J v = Vi -f 002'- — floP, = y v — Xx t 



VORTEX CONCEPTION OF ROTOR AND MUTUAL 
EFFECT OF SCREW/PROPELLERS 

A.M. Lepilkin 

In existing studies on the vortex conception of a screw, a /77* 

screw system with an infinitely large number of blades is used (which 
permits avoidance of the extremely complicated accounting for the ef- 
fect of nonperpendicular running off of the vortices from a rotor 
line and the time variable circulation of airspeed around the blade 
cross section), and only steady state conditions are considered. 

G.I. Maykapar in 19^7 [1] and A.P. Proskuryakov in 1956 [2] 
used a system with constant circulation over time. G.I. Maykapar 
generalized the method of Zhukovskiy for the case of a slanting vor- 
tex cylinder, and he proposed breaking down the vortices into circular 
(parallel to the blade planes) and rectilinear (along the cylinder 
generatrices), and he gave formulas for inductive velocities only in 
the planes of the blades in the form of definite integrals. A.P. 
Proskuryakov considered only a flat vortex system and broke the vor- 
tices down into longitudinal and transverse, and he gave formulas for 
the axial inductive velocities involving elliptical integrals of the 
third kind, which changed in complex planes, the parameters of which 
are determined by the roots of an equation of the fourth degree. 

A vortex theory of the screw is proposed below with a variable 
(according to the blade and by its azimuth) circulation. A three di- 
mensional vortex system is used as the base, the form of which takes 
account of the conical nature of the relative locations of the blades 
and the first harmonic of their flywheel motion. The results obtained 
are suitable for all normal conditions (with the exception of the 
case of axial movement of the screw against the force of gravity). 

In Section 1, formulation of the problem is given, and expansion 
of the circulation in a Fourier series by blade azimuth is proposed. 

*Numbers in the margin indicate pagination in the foreign text. 
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In Section 2 and Section 3, equations are given for the inductive 
velocities in space, in which in Section 3, it is for a screw with an 
infinitely large number of blades. 

In Section expansion of the inductive velocity of a screw by 
the second blade azimuth is given. 

In Section 5, multiparameter, improper integrals are presented 
in the form of a combination of elliptical integrals and elementary 
functions . 

In Section 6 (Appendix), a method is shown of reduction of el- 
liptical integrals of the third kind with a complex parameter to in- 
tegrals with a real parameter. 

1. Formulation of Problem 


In steady state motion of a screw, angle e between the blade /]_S 

axis and the plane normal to the axis of rotation of the screw changes 
periodically 


8 = a o + a ie cos Tp -f a u sin \{> -{-••• ( 1 ) 

The axes of the blades (if the short distance of the hinges from 

the screw axis and the higher harmonics of flywheel motion of the 

blades are disregarded) form a circular cone with angle of taper a Q , 

the axis of which is deflected from the screw axis by small angles a-^ c 

(in the plane of movement of the screw axis) and a, (to the right or 

1 s 

to the left ) . 

We introduce a clockwise coordinate system with origin 0 in the 
center of the hub of the screw by drawing the y axis upward along the 
axis of the "blade cone" and by directing the x axis forward so that 
flight speed vector V is in the xy plane (Fig. 1). We define blade 
azimuth ip as the angle between the projections of its axis on the xz 
plane and the x axis. In calculation of the forces and moments de- 
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veloped by the screw, it will be easy to change to a similar coordi- 
nate system with the Oy^ axis along the axis of rotation of the screw. 

Angles ag, a^ and a^ g can be calculated 
approximately beforehand based on axial induc- 
tive velocity v^ averaged over the propeller 
"disk" (it can be determined from for example 
the theory of Glauert [ 3 ] ) - 

Because of the smallness of angles a^ and 
a ls , the angular velocity of rotation of the 
"blade cone" can be considered uniform and 
Fig. 1. equal to velocity Q of rotation of the screw 

around its axis. In analysis of the mutual ef- 
fect of the screws, the axes of the "blade cones" can be considered 
parallel. The blades of the screw are assumed to be uniformly ar- 
ranged and identical in all respects. 

In determination of the shape of the vortices flowing away, by 
disregarding the projections of inductive velocity v on the x and y 
axes, which are extremely small compared to the velocities of points 
of the blade V and Pr, each element of a vortex can be considered to 
move parallel to the y axis from the point where it left the blade. 
In order for such a (linearized in shape) theory to retain meaning 
up to hovering (V x =0) inclusive, velocity V^*=V y +v d of departure of 
the vortices from the xz plane averaged over the "disk" of the screw 
(or in the xz plane) should be introduced. 

Steady state conditions, circulation r of the airspeed around 
the blade cross section at radius r depends periodically on blade 
azimuth \p . 

We introduce the concept of an abstract screw of radius p with 
constant circulation y(<J/) along the blade, and we determine corre- 
sponding inductive. velocity w. 
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Let t be the interval of time from a given moment to the time a 
vortex element flows away from the blade. The blade vortex connected 
with circulation y(i p) and two "end" vortices with circulation y,= 
y(ip-nt), flowing away at radii p and 0, form a closed circuit filled 
(because of the law of conservation of circulation) with "transverse" 
rectilinear vortices with circulation dy = (dy./d<}> )d<j> , the ends of 

<p $ 

each of which rest on points of the end vortices which correspond to 
equal values of t. The coordinates of the vortex element in the 
x,y,z system will be (Fig. 2) 

l = p' cos(t|>— 9 ) — 9 -^ -n = a 0 p' — ^ qj ( 11 ) 

2 = — P' Sin — 9 ) (<p = n/) 



Fig. 2. 


Only <}> changes along the out- 
side (p’=p) and inside (p*=0) ends 
of the vortices (0<4><°°) and only p T 
changes along the transverse vor- 
tices (O^p’^p). The equation of the 
connected vortex will be <|>=0. 

We call the p ,= const line 

(which is transformed at V =0 into 

x 

a screw with pitch and is pro- 

jected on the xz plane in the form of 
a trochoid at V^O) a helicotrochoid. 
The p’=0 line will be rectilinear to 
the axis of the vortex system. 


The inductive velocities which correspond to a given vortex line 
with circulation f are determined by the integrals along this line 

r * = 4n^ r (a^T d S- i-r d v) 

i^W'sr'r'' 5 - 

A- £5 3) 

/= V (*— 6) l +(y — nP + f* — 5)> 


( 1 . 2 ) 


where £ , n and £ are the coordinates of element ds of the vortex line. 
These integrals can also be applied formally to vortices with circula- 
tion of variable length with such integrals considered to be branching 
"transverse" vortices. For construction of the velocity field in the 
xyz coordinate system semifixed to the propeller, Eq. (1.2) can be 
utilized, since their result depends only on the differences of the 
coordinates of the points of this system and the points where the 
vortex element is located at a given time. 

At a fixed point of the xyz system, the inductive velocity will 
be a periodic function of blade azimuth ip with period T=2ir/k, where 
k is the number of blades. 

A change to a multiblade screw can be accomplished by sum- 
mation of the inductive velocities of a single blade screw shifted 
by ip in t, 2t, . . . Therefore, only the cases k=l and °° are con- 

sidered further. 

For a change to a screw of radius R with variable circulation 
T(r,^) over radius r, the following must be used in equations for the 
inductive velocities of an abstract screw with circulations y(^) , 

T (*> = - {£},_ * d-3) 

and the integral must be taken along the blade axis (0,<p^R). 

In order to simplify practical application of the theory, the 
following expansion should be used 

oo 

T ( r » S’) = F 0 (r) -R ^ (Fnc (r) cos m{> -f- r n » ( r ) s * n ™J>] ( 1 . l \ ) 

n=i 

A screw with unit complex circulation 

Y n =' e 1 "* = cos nip + i sin mf (n = 0, 1, 2, . . .) (1.5) 

will correspond to complex inductive velocity 
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w =w +iw 
n nc ns 


(n=0,l,2, . . .). 


( 1 . 6 ) 


The Inductive velocities of circulation harmonics r o (r), 
r nc ( r ) and r ns (r) w111 be 


/80 


V ”° (x, y, z) j Vnt (Xy y, Z ) = _ J tt ^ ( p)^J. dp |(l . 7 ) 


2. Inductive Velocities of k^°° Screw 

By applying Eq. (1.2) and (1.1), passage around the edge of the 
vortex surface must be carried out so that the connected vortex 
passes from p’=0 to p'=p and that the "transverse" vortices have to 
pass in the same direction. The following equations are obtained in 
this way: 


P 

« [sin (t - f) £ | + «0 4; f ]d<, dp’ - r (« j (sin * | £ + «,„ £ ') i? / 

o 0 

co 

4^= - P J. T<p [cos(t|)-9)±i_ s in( 1 l J _q,)A|] d( p + (2.2) 

0 J . 

P CO 

+ S S w [ sln -k T + cos (* - <p) £ T 1 dp' - 

no' -J 


- TT S T. j? (4 - -i-) it + T M>) S (si. * i -i- + cos ± ±) dp- 


4n». - - P $ r,sm( + — <p) A i + £ J r. £ (| _ ^ _ 


(2-3) 


0 0 


V* 


M r . Mi; - -c)^— 

0 0 

** = [x — p cos (^ — <p) + <pV x / Q ] 2 -j- (z + P sin — 9 )]* + 

+ (y — a 0 P + #v* / Q) s 


(2.4) 
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Here , Z * and Z are the values of £ at p=p* and p=0 respective- 
ly, and £ q ’ is the value of l' at cf> = 0. Differentiation with respect 
to x, y and z is performed before integration with respect to <f> (the 
condition of convergence of improper integrals). It is easy to see 
that 


i!W " = P 5 r, (|~ T + i I) d t ~ Tt f r, 4- (t - T*)< *p + 


0 0 n 


(2.5) 


By substituting Y=e ln ^ and Y (J) =e in ^ < * > \ complex functions of 
the type of equation (1.6) can be obtained for each projection of the 
inductive velocity. 

3 ■ k=°° Rotor 


As k-*°° and with preservation of quantity kT , period T=2ir/k of in- /8l 
ductive velocity v(x,y,z,t|;) with respect to ip tends toward zero. In 
the limit k=°°, function v does not depend on blade azimuth and evi- 
dently equals the average value of the inductive velocity of a single 
blade screw (with the same circulation) in one rotation of it 

Tt 

v° ( X , y, z) = ^ V ( x , y , z; (3-D 


The use of function kv° for a k^°° screw designates approximate 
solution of the problem, the inaccuracy of which decreases with in- 
crease in number of blades. A k=°° system, as numerous experiments 
have shown, gives satisfactory results with k^3 and small values of 
V/J2R, when the coils of the vortex lines are located quite close to 

each other. Rotor parameters V /flR and V */ftR also are small (usual- 

x y 

ly less than 0.3). 

By applying a' k=°° system in determination of circulation r(r,i|>), 
a blade can replace a vortex line since, with k=°°, function w can be 
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integrated everywhere with respect to p. 


With k=«, the vortex system of an abstract screw consists of a 
vortex cone (connected vortices of the blades) and a vortex shell, 
each cross section of which is a circle in the y=const plane. The 
space inside the shell is filled with "transverse" vortices. In each 
infinitely thin parallel xz plane, the circulation distribution layer 
of the shell around the circumference is identical, if the calcula- 
tion is carried out from the same helicotrochoid ^=const. 


In determination of an Eq. (3-1) type function for each projec— 
k-^-on of inductive velocity w, a detailed accounting will only be given 
for Wy. Transformation of the double integrals should be performed 
ahead of time for y(^)^const . After integration over (j> by parts, 
there will be 


oo • 



(3.2) 


In determination of the average value of Eq. (3.1) for a period, 
variable a =tp-<p can be introduced for each fixed value of <p . Then, 
a=ip can be assumed. The hydrodynamic meaning of such a transforma- 
tion is that summation of the action of the vortex elements in each 
parallel xz plane of the thin layer is carried out from the same 
helicotrochoidal surface <|;=const. 

In this way, we obtain 


o — * 0 —It 

o 00 -« 

Z* = (x - P cos t|> + 9^)*+ (Z + pain + (y - a 0 p + <p 


( 3 - 3 ) 

( 3 . 4 ) 
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Use of the equalities 


a j_ n a 1 V v * a 1 

dx l ~ V x d<p l V, dy l » 


&T**‘ 


lo 


(the first is necessary only for the w transformation) gives 


Aitw x ° = ns - S ^ S ^ (t ~x) dlt,<f<p + p \ l Sr C0S '* ,dl * ,d< P ~ 

O — n 


~ 5 fir l &' 003 ♦ 1 d *' d< f w \ iS' ! sin * 1 tf ♦ 1 rfp ' ~ a ° \ir \ sS ^ rfp ' 

0 0 — * 0 — R o — ji 

* — n. 0 — R 


a * 

P n 


0 — r 


0 — R 

7 . 


1 ( 3 . 5 ) 


(3-6) 


0 — R 


0 — R 

+tj 5 l ******<?+ 


(3-7) 


— n 0 — r 

+ Sf|r5 l ygcostj>^rfp' + a 0 ^^ 

o o — R 0 — R * — * 

(/.*=*» + »• + **) 


0 — R 
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tions 


There now should be introduced variable £=x+<f)V /ft and nota- 

X 


X = V v * / V x , y, = y — a 0 p — Xx 
P = to* — 2© (5 cos rj> — z sin ^) + l 2 + z* + (y» + *•£)* 


Then , 


a i _a_j__ /_a_n 
a P / +4 °ay i — la«o i j u _ P 

4™,° = - X [K t (p) - K z (0)] + ^-AT V (1) ( P ) -^^ v (1> ( P ') dp' - 

0 

-^ 5 ^ sin ' i, ^ dp, “ ao i'£' Sz^wp' 


(3.8) 

(3.9) 


( 3 . 10 ) 


0 — R 


0 — R 


(3.11) 
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( 3 . 12 ) 


— v ? x K* (P) “ K* (P) + Kt ( 0 ) + ^ "Jr dtydp' 

0 —n 

4nw z ° = (1 + X*) [K v ( P ) - K v (0) ] + X J 2^ (1 - JL) _ 

— n 

-yZ K v (1) (P) + (P')dp ' $ g^cos r}> *|>dp' +«o $ A J m d y tfp'j 


0 — n 

*•*- '*■+»•+** 


Here , 


| (3-14) 

X — rt 

.V-S*!®**** 1 (3 - 15) 

X — 7i X — n 

l* = p* — 2p (x cos i{> — z sin tp) + x* 4- z* + (y — a 0 p)* 


(3.16) 
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Equations (3-14) can be presented as 

*„«=- j *<*♦. *.«■> = - j dgd* 


(3.17) 




p — ^ cos ^ + z sin tj) 


d£drp 


Functions K ^(p) and k/ 2 ^( p ) are obtained from K (p) by 

y y y 

replacement of y(^) by y(ift)cosijj and y (^)sinij;. Since 

P = A + 2BI + C£* 

= (z + p sin ij)) 2 + P* cos* ij) + y # *, B = Xy e — p cos \|>, C = 1 4- X* (3*18) 


then 


f d\ _ V C B + xC r | d\ _ A + xB B 

) I s ~ H HD 1 ) !• “ 


77n =W , H = AC-B* 

HD HVC (3 . 19) 


H = (z + p sin tp) 2 (1 + X 2 ) + ( y t + Xp cos ip) 2 
D* = p 2 — 2p (x cos tp — z sintp) 4- z 2 4- z 2 4- (y — a 0 P) 2 


Consequently , 


KM - - \ ^ -.<■*> + (^ - 4 ^ <3 -20 ) 
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&(?) = — ^ lM(/C-^tf£)(z + psinr}))£i^ 

— « 

ff«(p) = - j {(Kc-^^) (P + 2 sin 'JO — [~yy— — :pr=) cos f}<ty ; (3 * 20) 


The remaining integrals with respect to ip will be 

* n I 

C 1 if) _g_ _j_ Ah f TW *cosiJ> 

J 2 n 0 * /o' a V~ ) £5 

K — n 

JL [ Ilf) W,i, _ [ T (f) * — P'cos tp , 

dx ) 2nl a - a V - ) -2JT 5* 


n 

$ 


. lif) Will = f I 'V) ^Tfamy > . 

dz } 2nl 0 ' ai * > ) 2n ® t 


T OW z + p'sin if) 
/o' 3 


i $ is&rw* 


( 3 . 21 ) 


All the integrals with fixed function y(ip) can be found nu- /8U 

merically . Calculation of integrals (3*20) is complicated however 
by the extreme peaks of the subintegral functions at small values of 
H. The case H=0 and D=0 is possible with y — a 0 p = y r p r — z i . 


The calculations are substantially simplified if the circula- 
tion is expanded in a Fourier series with respect to \p. For 

Yn= ein ^ (n=0 , 1, 2, . . .), integrals (3-14) are transformed into 
special functions 


*»>(?) = 


x — :: 

^n+l. y(P) 3 * K n-1 . v (P^ v (2) 

2 ’ ^ nv (P/ 


(q = y, z, CD) 

*n + l.„(P)-*„-,.„(P) 
2 i 


( 3 . 22 ) 

(3.23) 


The substitution il»=ot sgn z should be introduced here. There 
will then be 


«•-» .1*1 


(0) 


x — rs 


(3.24) 


P = to* — 2co (£ cos a — | z | sin a) + ,g 4 + z a -f- ( y e + A.£) a 


11 


Functions of the type 


where 


F » **’ g ' h; \ \ lni dad ^ (“>o. *>o,a>o), (n = i, 2 ... .) 1(3-25) 

X — « [ 

i 

! 

^ — 2to cos a — g sin a) + £ 3 -f g* -f. (h -f Jt£)* 1(3*26) 


and its derivatives over g, h and to (y=0, 1, 2, . . .) are discussed 
below (in Section 5). These functions are presented in the form of 
a combination of elliptical integrals and elementary functions, in 
which 

F * = F »*. + iF »" F m = F w + iF^ q (3*27) 


It is now evident that, with X g = | X | sgn (Xy ), 

K nq (to) = Fncq (x, I \z |, I y e |; to, \) + i sgnz F nsq (x, | z |, |yj; to, X e ) (3.28) 

Consequently, in the equations of Section 5, it is sufficient to 
assume g= | z | , h= | y g | , where | h+Xx | = | y-a^p | . It is then necessary to 
assume <o=p. 


For y=e 
integrals . 


in ^, the remaining integrals are reduced to elliptical 
Actually , 


3 f COS ... T„0M ,, . r e ‘ (n+l) * , i (n-l) 4> 

a f T„0|>) , P e 1 "* JU e Mn + i)*_ i(n-l)4, 


4n/ 0 /s 


cty i 


d f TnW) e ‘ iB * ,, , / f e ‘<n+l)«» + l(n-l)* 

to 1 2 ^ = — x \ 5 ^,^+P 5 ^ 

—7i — it — n 

f T„0M 3 1 . f d 7 „ JiJ) . P e in * , 

J 2n t>«y /o' ^ ~ J dip 2n/ 0 ' ~ ~ in ) 2 n/7 ^ 


(3*29) 


With the structure of expression (3*16) taken into account, it /85 
should be assumed that 


cos0 = x / r, sin 0 = — z / r, r = J/’x* + 


(3*30) 


12 


After this, use of the substitution \p-& = TT-2<p gives 

2 (-<)•«.<-«(*) ,, 4 P r 

J <’ » ((P + /■)’+ <J- (P + 

n/2 


P (-T71/2) 



cos 2n<p <fq> 

(i — k* sin* q>) m/2 



'(3.3D 

(3-32) 


For an abstract screw, integration must also be performed with 
respect to p* from 0 to p. The way of calculation of the correspond- 
ing double integrals can be indicated. For a real screw with vari- 
able circulation however, r(0,ifO = r(R,i|O = 0 . Therefore, integration 
over p* proves to be superfluous. Actually, for a g(p) type func- 
tion with r(o)=r(R)=o, 

R P R 

J (P') d P') d P = — (P) r (P). d? 


Function 


a 

E n * v> (At, a) = ^ (1 — k 2 sin 2 <p) v cos 2n<p.rfcp 
o 


(3-33) 


can be called a generalized elliptical integral, since 


£o ( - ,/,) {k, a) = F (k, a), Eo { ' U) (k, a) = E (k, a) ( 3 . 3 *0 


It is easy to verify by differentiation with respect to a that 

I(n - V - 1) M V) (*. «) + (» + V + 1) E n V (k, a)] 4. (3 35) 

+ 2n (2 — fc 2 ) E n (v) (k, a) = 2 sin 2na (1 — A 3 sin 3 a) v+1 

The base here will be functions E Q ^ v ^(k,a) and 

E 1 (V) ( k , a) = - E 0 (,) (k, a) * JL/? 0 < V+1 > ( k, a ) (3-36) 


it also is easy to determine that . 

E 0 <v+V (*. a) - (2 - k 3 ) E t < v) (k, a) + 

•f- — 1-r (1 — **) £o <v ' 1) (*» a) = ** --77 (1 — ** sin* a) v sin a cos a (3-37) 

V ~r l Z (V 1) 
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where elliptical integrals of the first and second kinds (3-3*0 will 
be the base. 

2 

With small k , when recurrent Eq. (3*35) is inconvenient because 
of the effect of small differences (increasing with increase of n), 
it is better to use the expansion 

: (3.38) 

a 

(sin* (p ) 1 cos2n<fd<f i ( 3 - 39 ) 

0 

By using the corresponding indeterminate integral [*J, p. 155], 
we have 

[ (sin 1 <p)* cos 2 n<p dq> = (sjn 1 a )‘ sin 2 r»a _ _n_ f (sin , 9) « sin ( 2 n — 1 ) y d<f ; 

J 2 (n 4 $) n 4 *sj sni(p 

But 0 0 

sin /cep __ ( 1 4 2 [cos 2 <p 4 cos 4q> 4 • • • 4 cos ( k — 1 ) cp] (fc = 3, 5, . . •) j 

sin (p \2 [cos 9 4 cos 3<p 4 . . . 4* cos (k — 1 ) <p] (k =» 2, 4 . . .) 


On this basis, the recurrent formula 


a> = 

n, $ 


_ (sin 1 a)* sin 2na 


2/i 


2 (n -f ,) n + s (a > n + i (^l. • + 

a 

f. (a) = ^ (sin* <p)* dq> 


' ( 3 - * 10 ) 


is obtained, where f (a) is a known function which can be determined 

s 

by the recurrent formula 

2 tf, — (2s — 1) — (sin a) 2 * -1 cos a, /„ = a, 2/i = a — sin a cos a 


The substitution of /l-k^ tan <j> tan 0 = 1 gives 

= - ft ~ . 1 — A* sin 1 <p = 1 ~ ** h, v 

1 — A* sm* 0 1 — A* sin* 0 ( 3 • 1 ) 

n /2 ^ 

\ (1 — A* sin* (pf dq> = (1 — A*) v+Vl \ — 

J J (1 — k 3 sin 2 0) v+1 


where a and g are connected by the relationship tan a tan g/l-k^=l . 



For a=Tr/2 , integrals (3.33) will be "full," 

7T/2 

E n^ (*) = ^ (1 — A* Sin 1 <p) v cos 2n<p dtp / (3-42) 

6 ' 


According to the identity 

(sin« <P)‘ = [2 ^ D-* cos 2(s~g)+ ( 2 * )] ~ 

n /2 

O n .,= j (sia*«p)*cos2«(prf<p=(_l)n^.( j 2*J^ i o n ... = 0 (,<„) 

Consequently, for n=l, 2, 3, . . ., 


E n (v) (k) = 


r(^r{a)-(„; 1 )rt 2 )f + t; 2 )( 2 “ 2 +4 )(^r-f (3 - 43) 


For v<0, all the coefficients of this series have the same sign 


(-l) n , so functions E ^(k) do not have roots. 


For a=0 when 6 =tt/ 2, Eq. (3-41) gives 

n /2 71/2 

^ (1 — A 2 sin* <p) v dtp = ( 1 — A 2 ) v+,/ * . 




(1 — A 2 sin* q>) v+1 


(3.^) 


For v=y-l/2, it follows from this that 


(k) = (i - kr^-™ (k) (n = 1, 2 , 3, . . .) 

_ E o (Vt> (k) 

(1 — A 2 ) 2 ’ 


Eo { ~' U) {k) = , 

' 1 — A 2 ’ 


Eo <-,/,) (At) = 


(3-45) 


It follows from Eq. (3-35) that 


eJ-' m (k) = — 4n 2 - ** 


(k) - 2 £=± (A) (3.46) 


2» + 1 A 2 ” ’ ' 2n + 

Ey u) (A) = (A), A 2 £ , l <_,/,) (A) = 2ZT (A) - (2 — A 2 ) F (A) 

E, <-V.) m _ _ 4n 2 - A 2 


r “< (» 2T drr £•„<-*> (« - e„hu (*) ( 3 .'H7 ) 

(1 — A 2 ) Cj,* /,) (A) = ZT (A), A 2 (1 — A 2 ) E^~’ h) (A) = 2(1 — A 2 ) A’ (A) — (2 — A 2 ) £ (A) 


Equations (3-^6) and (3-^7) are necessary for calculation of /87 
integrals (3-32). 
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Point k=l is singular. As is known, 

^ ( * ) = ln Fr^ + T( ln ... I (3-48) 

It is therefore convenient to introduce functions limited to 
point k=l 

An (*) = (- i) n (l - a*) b™ (k), A n (l) = 1 ' (3.49) 

which are determined according to Eq. (3.47) by the recurrent formula 

A - 4n 2—A* . _ 2n+l . .. 

n+x 2n — 1 k * n 2n—i An ~ u ( A -n - A r) 

(3-50) 

A, (k) = E (At), k *A t (A) = (2 —*»)£(*) — 2 (i — A*) ^ (A) 


4. Inductive Velocity Harmonics and Mutual Effect of Screws 


We discuss a pair of screws with parallel "blade cone" axes, 
with origin 0 of the xyz coordinate system placed in the center of 
the hub of the screw, the effect of which is under consideration. 

We place origin 0 ^ of the identically oriented coordinate system in 
the center of the hub of the second (to be calculated) screw. As 
before, let ip be the blade azimuth of the first screw and 0 simi- 
larly be the azimuth to be calculated of the blade of the second 
screw (Pig. 3). Actually, in a two rotor helicopter system, a pair 
of counterrotating screws is used. This is then easy to take into 
account, by substituting the sign of azimuth 6 in the final equa- 
tions. Let L be the distance between the screw axes, 0 be the angle 
between normal L which connects the axis and the xy plane (Fig. 3) s 
the "angle of slide" of the pair. The coordinates of center Og of 
the hub of the second (to be calculated) screw will then be 

x 2 =L cos p, y 2 , z t = —L sin p (4.1) 

The coordinates of the point on the axis of the blade of the second 
screw at distance r from the axis of this screw at angle "of taper" 
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a Q2 of the location of this blade will be 

* = X % + r cos 0 , y = y 2 + a 0 2 r . z = z 2 — r sin 0 : (4.2) 



Fig. 3- 


The problem is formulated: to find 

for an abstract screw with circulation 
Y n =e ln ^ the expansion 

oo 

Wn (r,Q) = w n ,o (r) + ^ [wn.-mc (r) cos m0 + w n ,m» (0 sin mQ ] 

m=i 

(4.3) 

With complex harmonics 

®n,mc 1 Wnc.mc *f" 4w n j, mci W n , ms — Wnc t ms + IWrx, ,mi 

(4.4) 


for unit circulation Y n =eln ^ and using expansion (1.4) of circula- 
tion r(p,ip), the Fourier series 

CO 

(r, 0 ) = v„.o (r) + 2 [Vn.mc (r) cos mQ + o n .m. ( r ) sin mQ ] ( 4 - 5 ) 

i 


/88 


can be plotted for the inductive velocity of a real screw, by cal- 
culating the required number of harmonics by the equations 


(• dT 

V nc ,0 = — \ ttlnc.O ( r » P) dp. 


dp 


v 

nc,mc 

v nc t ms 



W 

nc.mc 

w 

nc,ms 





v 

ru.mc 

V 

ns, ms 



dT„ m 

tvs, me ns 


d P 


d? 


(4.6) 


for 0<r<R2* Usually, R^R (the screw radii are the same). 

In the case x^ = y 2~ z 2 = ^ 3 the harmonics of the natural inductive 
velocity of the screw are obtained. Subsequently, only an analysis 
of axial inductive velocity w y ° of a k=°° screw will be given (w x ° 
and w 0 are unimportant for the operation of rotors). Exponent 0 is 
dropped here. 

Equation (3.6.) for Y=e in ^ can be presented as 
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S£i*+- 

0 — n 

+ F 1 3r ( P ( , $ r ~ + a °i T-) ~ IT 4- (t - -^)}^ rf< P 


( 4 . 7 ) 


2 

where H function ( 3 • ^ ) a in accordance with Eq. (4.2), will have 
this form: 

P = (** — P costj) + r CO30 + (fV x / Q) 2 -f (y v + q>V tf * / Q) 2 + 

+ (z 2 + P sin tj> — r sinG) 2 , y v = y 2 + a<Yi r — a 0 o (4.8) 


Differentiation is carried out in Eq. (4.7) only up to the in- 
tegration operator; if n=l, 2, 3, • • the integrals for p=0 are 

absent (this subsequently becomes less evident). On the assumption 
that 

a 2 = P 2 — 2pr cos(>f — 0) + r 2 , e ix - p ~ re<(IMl) 


equality (4.8) can be transformed thus: 

P = (O' - 2 a D cos<>|> + Z - «) + Z>’ + ( Vv + 9 F„*/ Q)< 

D ' " (l * + ^ ' 0,1 + **• e " - <*. - / 2) / D ( 4 . 9 ) 


where 


It is significant that w and e 1 ^ are functions of i|<-0. In Eq. 
(4.7), the average value of the subintegral function for the period 
is used. It can therefore be written that 


z>in6 


to„ = -»«- J 4£ !! *i’7 1 + 

0 — 7t-0 

f T d 1 . v x a , i i .. 


OO 71—0 


and ip-e can be replaced by ip ; the result will be 


+ ~\ 1 £{> ( vt + * U) - £ 4- (4 - £)}<♦* 

Z 2 = (o 2 _ 2 co£» cos (t|> + 0 + X - t) + £>* + (y v + / Q) 2 
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6> 


- y p 2 — 2pr cos t|> -fr 2 , e ix = ^ — — 




(4.10) 

(4.11) 

(4.12) 
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We call this transformation the w transformation, 
the series 


Jn9 


do 


—j~ — d 0 + 2 ( a mc cos mQ + dms sin mQ) 

m=i 

f jne a r 


i(n+m)0 , *(n-m)9 

± = do 


2nl Y± 1 


We introduce 


(4.13) 


By introducing variable 0=0-a, where a=x-x-9 does not depend on 
0, it is easy to obtain 


P = to 4 — 2 «>D cos 0 + Z > 2 + (y v + q*V / Q)* 


(4 . li|) 


tv r 

sr<»> 


mi 


v* 


n)d - 

f* J 


<Z J(n+m)<3 « J(n+m)9 i(n-m)o /* i(n-m)8 

f Ssr-"±W$- 1 -'-‘» 

— Tt — K 


2nl 


or 


a 0 = ^o e ~ ir " 1 ', ' = e * 

a ™ K±1 


where, if the integration variables are selected from the left. 


71 

A ±m (p, (p; tp) == e *(n±m)(T-z)Tim^ ^ 


J(n±m)9 


2nl 


dO 


(4.15) 


Now, based on Eq. (4.10), it follows that 

P « ^ CO K 

^"0 = — *« \ \ A 0 (P't 0) — if ^ 5 "ST ^ P ’ ^ _ A ° ^ X 

0 — It o —« 

oo « 

(; x dx P d< P + P S 5 [l£ A ° (p. 9 ) + °o-^-A 0 (p, 9 )] <ty<fip («=* 0 ) 

!;, o — n 

■J c?i|) rfcp — 


Mj( 


3 ^ m (P.<P)±^_ m (P.<P) , ^ /l m (P.<p)±/1_ m (p.<p) 

op VTr +a °^ 7±r 


r r rd 

A m (P. 9 ) ± A -m (P. 9 ) 

d A m (°> 9 ) ± A -m ( 0 « 9 )"j 

J J 0 z 

0 — K u 

VTT 

oz V±T J 


After substitution of the order of integration over \p and <j> by /90 
introducing the notations 
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(4.16) 


V r p -in© 

= c<|lT \in d & d< ? (±^=1.2.3....) 

0 — it 

y CO ” 

K m (a) = A jeV* ^ d<p (± n = 0, 1, 2 ,...;? = y, *, <d) 


since D=L and x=B for <J>=0 we obtain: 


for m=0, 


r r r* nX# r 

4^.0 = - ine ^\ l -2*7- S S*v <W *Mp' + : 

0 — n —it | 

+ fr { S + “« S-Tir^w^H <4 - 17) 

— n — n / 

~ ) + S ± r Kn * {r)d ^ 


for m=l , 2, 3, . 


4n 


“'n.me . f r P 

“’n.m, m l jATTi J J 


«-i(n+m)X'-im4> p l{n+m)& 


0 —« 


2 Jtp' 


) 2«V ^^ dp,± 


± TWr 1 1 — i,p- $ TaT-w^M + 

0 —it — n 

it 

Qp r p tf -Mn+m)x-im<J» r a „ _ ! 

+ v^{ ) 2n 7Tr'F n+m ' w ^ ~ 1 + m ) ^ n+m (®) ^ i 

—it 

« 

p c -t(n-m)X+im4» r ^ 

± _i ~ 2n K TT r— (®) i (» - !») *n- m (CO) g]*|>} + ( 4 . l 8 ) 

* 

n n 

, Qp f p c -t(n+m)X-im<J> « _-i(n-m)X+im^ . 

+ °° VZ l J 2,-t ^TT Kn+m - v (<0) ^ ^ ' 2k V^tT Kn ~ m ‘ v (&) ^} “ 

71 — n r 

n . _ 

f r tf -»(«+Tn)X-im^ p -i(n-m)X+im>V , 

~ l J 2* Ffr * n+m>1 (t0) ^ ± \ "TnVTT (co) M + 

* — 7t r 

U " ^- i ( n +™) x a -i m 4' ” -i(n-m)X a +imcV 

_ 2n V^f Kn+m - z (r) ^ * S e 2 «VTr (r) M 


Here and subsequently, x* and w 1 designate quantities x(p) 
of Eq . (4.12). with p=p' 
follows from Eq. (4.12) that 


w(p) of Eq. (4.12). with p=p', and x is the value of x with p=0. 

cl 


and 

It 
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cos 


P r cos tf ) __ rsin^ 5co 

* CO * ^ (0 * dp 


C0S X» 


*x. 

dp 


= - S -^ (4.19) 


where x(~^) =_ x(^)* Therefore, in the notations 


V _!£t ± A* 

~ da* a A *» 


r + aA '° - K - 
i^o = ^j--A 0 


(4.20) 


the following relationships are obtained: 


for m=0 


da dr Jx ,-ix 

^. sr -,>/f, 5 i=4 r v+ W 


2 ^ 


4 “ v ” = - $$ if* 5 ^ <» ^4/ + 

0 0 — n 

+ 77{tS [cos < n ~ Dx^n + ((o) + cos (n + 1) X ^(o)] d* + \ 

• - 

+ a 0 ^ cos n%K nv ((xi ) ^ cos n%K nz { <•>) + ^cos n XaK n z (r) <fi{> 


(4.21) 


for m=l , 2 , 3 » • • . , 


4. eiin ~ mYfi C f cos [(n -m)x'- m* I f e 1(n - m)ft jajuj /| , 

± ^T)) 7 ) -2^<K>*|uZp} + 


s- 4 n a n - TO ^ — ,-n J l__ Cf cos [(n + m) x' + m^] r e f < n+m >* 

\» iy-±r jj p' i ~2^r Md * dp ^ 

e Un-myj 3 

"yTr 

V V, — Jt j 

I 

k 2T^fi cos( (' 1 + “- *)x + »"H Y§r i * + \n.22) 

{ * '_■ 

+ 5 cos l( n + m + !) X + ™I>1 d y ± 

0 * 
n 

\ cos l(n — m — 1) x - w»|>) ± 

o v * * i 

rt . . __ 

V ± \cos [(n — m + 1) x — m>Jj] -f 

o 1 : 

? n 

p" ° 0 ^ cos + m ) X H- wijj) ± 
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db^cos [(n — m) x — mr|>] - j: y j ==- — 

0 

cos [(n + m) X + m\j)] cty ± 

± ^cos[(n — m)% — m^] -2^==^-^} + 
0 “ 

||j cos [(n + m) xa + mrj)] ^r= ~ ^ ± 

cos [(n — m) Xa — w'i’l ~ ‘ty} 


It should be remembered that, for n=l, 2, 3, • • * , the in- 
tegrals for p=0 (o)=r) are absent. A change in Eq. (4.16) to the in- 
tegration variable $=x~+<J>V /ft with notations 

d X 

X = V y* / V xt = y? Xx 2 


C 1 * .23) 


in accordance with Eq. (4.14) and (4.9) 3 gives 

d'Q 


r p sv-oav 

K „(©) = \ 4- \ 7 t= = rdl (y=f 0) 

J 2n J Y(^ — 2aD cos Q + D* + (»/,+ XS)» 


X, —It 


D = V¥~+ z 2 2 . e" = (5 — *2j) / D 


(4.24) 
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By introducing the substitution 0=a sgn Z 2 3 it can be determined 

that s 

^ (w) = p^exp(^^ da£iS 

P = to 2 — 2© (l cos a — | z 2 1 sin a) + £ 2 + z 2 2 + (y m -f X£) 2 ( ^ • 2 5 ) 

Thus, Eq . (3*25) introduced above and functions F^(x,g,h; to , X ) , 
(y^O) described in Section 5 arise. With X#=jX| sgn (Xy # ), 

K v .{(a) = F| ll | C (x J ,|z 2 | t |y*|;<o, X*) + isgn (ysgn z a ) (x 2 , | z 2 |, |y # |;© , X*) 

(±H = 1,2....) (4.26) 

AT m (o) ) = /Wa(x 2 ,| z i|. I y*|; ©.X*) + i sgn (y sgn aQ (ar a ,| z 2 1, | | y+ 1; ©, X*) 

(±y = 0. l. ...) 

where q= | Z 2 | 5 y* > to . It is therefore sufficient to assume in the 
equations of Section 5 that 


22 


x = x % , g — |z*|» h — |y*|, | h + %x\ — |y v | 


Equations (4.21-4.22) also include the integrals 



— n 


\ 


£ V-*d§ 

2 n — 2 coD 0 C03 d + DJ + y v 2 ’ 


jD 0 = L G 1 — 0* 1* • • •) ^ ij ^ 2 7 ) 


By operating as specified above in Section 3, the following can 
be obtained: 


[ (~D* 

2n/ 0 - * y^-\-L)' + y v * 


E^(k), 


k 2 = 


4o)L 

(o) + ^) a + y^ 



cos 2 fi 9 d(\ p 
|^1 — A 2 sin* (p 


(4.28) 

(4.29) 


Functions ^ 1//2 ^(k) were introduced and studied in Section 3. 
We study functions to and x- According to Eq. (4.12) and (4.19), 


CO 


j/"(P — r)* + 4 P r sin a -|- , X = arc tg - ^ 

co = p, X = 0 with r = 0 
co = 2r sin 7^, X = 7*(«— W with p = r 
co = r, x = n — r l J ' with P = 0 




(p>0 

(p = r) with 4* = 0 
(P<d 


(4.30) 

(4.31) 


It is easy to compile tables of functions w(t|>) and x('P) For a /93 
number of calculated radii, by taking the value of ^ every 10 or 15° 
but a smaller step with small values of ip. It is evident that 
0 <co<R+R 2 ; usually, R 2 =R. R=1 can be adopted in the calculations in 

compiling tables for O^r^l and 0$p^l. The author has tables of oj ( ^ ) , 
x(4>) and cos [(n+m)x+m^] for +m=0, 1, 2 and n=0, n=l. 

With r=0 or p=0, to does not depend on ip . Functions K^co) and 
K (to) can therefore be taken out of the integral sign. With r=0, we 

r* M. 

have x = 0, <o=p, and it turns that w =0, but 
A * n,m 5 
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(*». 32 ) 


r r , inS do' 

4 ™„.o ^ ^ T + 

o — * 

+ &* {Knu (p) + aoK nv (P)> - Knz (p) + Knz (0) V = °> 

* X 

In the case of p=0, all w so. 

5 n,m 

The practical value of the u transformation is that, with a 0 =0, 
parameters r and p are combined in one parameter to (4.12) which, for 
a given relative location of screw pairs (x 2 , y 2 , z 2 ), with a given 
orientation of it relative to the speed of movement (angles 8 and 
arctan X), is the unique argument of functions F (x, g, h; to, X). 

Even with ag^O, the first three parameters of functions F^(to) and 
F (to) provide complete solution of the problem, since it is suffi- 

r 4 

cient to assume h= | y* | = | y^-Xx 2 | . 

After calculating the required number of functions of to(t|0 for 
the lower harmonics, integration with respect to ip can be performed 
by numerical methods. 

The following individual cases are of great importance: 

1. z 2 =0 (flight of a "longitudinal" pair without slipping); 

2. x 2 =0=z 2 (pair of coaxial screws or screw with a Q ^0); 

3. x 2 =y s =z 2 =0 (isolated screw with a Q =0); 

4. x 2 =0=y* (flight of "transverse" pair without slipping). 

In the latter case, the functions have two branches each: 

(4a) X 2 (1 - 4 . X s ) — a>*> 0, (4 ) L* (1 -h Ji 2 ) — co*<0 

Case (4a) occurs in the absence of "overlap" of the screws 

(l>r+r 2 ). 
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In Section 5, functions F (x, g, h; w , X) are written out for 
(1) 'g = 0, (2) g = 0 = x, (3) g = * = h = 0, (4) x = 0 = h 

In the first and fourth cases, the parameters of integrals of 
the third kind are real. In the second and third cases, functions 
F (w) and their derivatives are simple algebraic functions. 

In the case g=0, the algebraic functions in F (w) and F (oj), 

y yq 

with Ah<0(Ay x <0), have two branches each with branching point h=o) | X | . 
This must be taken into account in integration with respect to if in 
Eq. (4.21), (4.22). Transition point if„ is defined according to Eq. 
(4.12) as 

• “ Ul = ly.l, 2pr cos if A = P 1 4- r 2 — (y. / X) 1 

in which, in segment [0, if„) , we have | p-r | y*/X | . 
that 

<■> > I y. / X|, 4> A = 0 with ( P- r)* > (y. / X)' 
u I y. / X |, if A = n with (P "F r)* < (y. / X) 1 


(4.33) 

It is evident 

(4.34) 


Transition point if^, in case (4), will be at io>L ("transverse" 

pair of screws with "overlap"). It is determined by the equality 
2 2 2 

L (1+A ) =oc . It is sufficient therefore in Eq. (4.33) and (4.34) to 
replace (y#/A) 2 by L 2 (l+A 2 ). 

In the second and third cases, it proved to be possible to cal- 
culate all the integrals of if by presenting them in the form of a 
combination of elliptical integrals and elementary functions. These 
results will be the subject of a separate article. In the fourth 
case, half of functions F^(w) and F (w ) are simple algebraic func- 
tions, namely: 

dF »c / dh and / dg with p, = I n ± m | = 1, 3, 5, . . 

F \lc> dFp, / Oh and OF^J dg with p, = | n | = 2, 4, . . dFjdg 

Integration over if also is successfully performed here in a 
similar way. 
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5* Special Functions of Vortex Theory 


We consider function (3-25) of five real variables 

°° n <iia 

F * (*» & h > *) = ^ \^[ dad % (n = o,i, 2 ,...) (5.1) 

x — n 1 

/* = (■)*— 2co (l cos a — g sin a) + I* + g* 4- (A 4- X?)*, to >0, * > 0 , A > <f( 5 • 2 ) 
After integration by parts over a, there will be 


7 t . 00 

A\lo. 


^ I . A n P _ + * c P sa dCt ' P=A + 2 Bt + Cl* 

—71 X 

A = (a* 4* 2co£ sin a + g 1 4* B = — to cos a, C = 1 4 ^ 


(5.3) 

(5.4) 


Consequently, Fq(co) does not exist. However, derivatives 

dF„ ' “ 


dh 


- — { rfS da, C g ± to sin a 

J 2jx J / dg- J 2 jx j (9 


da 


then 


d/V 


dF tt 


d I» = - 

* oo 

r e'» a r 

c»> — £ cos a 4^ g sin a 

dco 

J 2n ) 

— It X 

p 

< 

1. Since 




= V"C _ B +xC 

1 

ii-p 

w/> 

A + xB B 

H 

H HD ’ 

i p 

HD HYC ’ 

D* — co* — 2u 

(x cos a — 

g sin a) 4- x 1 4- g 1 

4 (A 

— 71 

B + xC 

\ „ ,M + *B 

D 

J g cus a y ^ D 


— 71 

n BA-xC\. . M + rfi 

B \ 

/> 

r-n d ~ 

Yc) 

-$&(*> 


•) (g *4 (o sin a) <fa 



(5-5) 


(5.6) 


Vc 






(5-7) 


— ^ ^ - ( e<tta f(vn !LLz £. ), , . , [A + xB B\ ) 

5(o J 2ntf IV Z> j (“ 4- g sin a) — ( — ^ — - — J cos a| 




II =(fh sin a)* (1 4 . X 1 ) 4- (A ^ X<o cos a)» > 0 (5.8) 

With oj= 0, only these functions are different from zero: 


r 


dF 0 
dg 
dF x = 

a* 

_a/l 
dco = 


— 


AJt 4~ xC\ 


o = g 2 (1 4 X 2 ) 4 A 2 , 

Do' = 

* s 4 j 2 4 (A 4 X*) 2 

/ g 2 4 A 2 4 xAA 

XA \1 

(5 

[ Do 

}Pcij 

l g 2 4 A 2 4 xAA 



V Do 

Kc/J 



For oo^O, function (5*8) should be presented as 


H = 


P (to cos a) 


(e — (O sin a) 2 (1 4. A 2 ) 4 (A 4 Aco cos a) 2 


( 5 . 10 ) 


In Eq. (5-10), /95 


P (x) = (x 2 — 2A (A 4. ig V i 4. A. 2 ) x-( H ig V~ 1 4 X 2 ) 2 — to 2 (1 4 X 1 )} X , . 

^ (** 2X (A iy V 1 4 X 2 ) x — (A — ig y 1 4 X 2 ) 2 — to 2 (1 4 X 2 )}, x = to cos a 


Equation P(k)= 0 is broken down into two. The first of them 

:< 2 -2 X(A4 igVT+tyx-Vi+ ig VT^J?)* — co*(i 4 a*) =0 (5.12) 

has the roots 

** = X (A 4 ig Y l 4 A 2 ) ± sgn A 1^1 4 A 2 T^co 2 4 (A 4 ig Y 1 4 A 2 ) 2 (5.13) 

and the second has roots and k^, which are complexly conjugated 
with roots and • It is easy to determine that 

1 _ R (<a C03 g) — 2 g (1 4 X 2 ) co sin g (5.1^0 

H P ( o) cos a) 

R (x) = A 2 4 2XAx - x 2 4 (g 2 4 “ 2 ) (i 4 X 2 ) (5-15) 

1 _ y -R (*<) — 2g (1 4 X 2 ) <0 sin a w cos a _ xi -R (*v) — 2g (1 4 X 2 ) to sin a 
1 H “Jj (to cos a — x v ) P' (x v ) ’ " (co cos a — x v ) P' (x*) Xv 

j G>_sinjj_ _ .4, P (x v ) to sin a — 2g (1 4 X 2 ) (to* — x y *) 

i H ~ 2j (to cos a — x v ) P' (x v ) ’ ' 

I V=1 


For and k^j with the use of equalities (5.12) and (5.13) s 

Z “!={(* * ** v T Tv> rr+xii [xi K«» 4 (a 4 ^ v;ftt 2 ) 2 )* 

is obtained and if, for these roots, there are introduced the radicals 
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- = igVl + X*)/l -PX*±UI Vco* -HA + .> /l + X‘)» 

V Xl* — co* 


( 5 . 16 ) 


R(i<)=-i2g(l+X 2 )/K 2 -w 2 . Therefore, by assuming that 

= (/= 7 ) v V V- o«, (/ZTT )v== / * for andK 2 i ( 5 - 17 ) 

'• _i for K,andK 

J 4 

it can result that R( i<)=-2g(l+X 2 ) /o 2 -k 2 and consequently. 


JL ^ r V* 6 * 1 x »* 4~ M sin a 2g (1 + X») 

// Zi v (co cos a — X v ) P ' (x v ) 9 1 * i 5 ' (x v ) 

W=1 


cos . ■ . Y T Vu* — x v * + to Sin a 

H sin ~ 2l v (tocos a- x v ) />' (x v ) 

v=l 


/©* — X V 1 


( 5 . 18 ) 


Paired equalities are used here and subsequently. By differ- 
entiating first polynomial (5-11) over k and substituting the values 
of and ic^ in P’(k) and utilizing equalities (5-12) and (5-13) > it 
can be determined that 

**'(*) = ± i sgn X 8g (1 + X*) Vu'-±(h-\r ig VTrW* (UX * l ) 

r v = - ( -l) v »-sgnX (v =• l,v =2) i (5.19) 

4 (A 4- Xx v ) K + (A + igVi-k A*)* 

and, for v=3 and v=4, it is sufficient to change the sign in front 
of i. 


Roots k can be complex. We will therefore consider separately 
v 

the real and imaginary parts of integrals (5-7), based on the nota- 
tions 

^ ( 5 . 20 ) 

In using identities (5-18) in Eq. (5-7), in place of A and B, /96 
it is sufficient to substitute the numbers 

A, = g* 4- 2gVu* — x v * -p (0* + A 1 , B v .= XA — x v (/ly C — By* = 0) ( 5 . 21 ) 
for which C(A v +xB v ')=B v (B v +xC) . After noting the identities 
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V to* — x* co sin a 
to cos a — x 

V to* — x 1 4~ to sin a 
o) cos a — x 


to cos 


n „.. f y<n»-M»-fr 


ci) sin a 
a) cos a — x 


•b V to* — x 1 + (0 si a a 


to sin 


a = Vrc* — X 2 4- (■> sin a 


o) cos a — x 


Co cos a — x 


and introducing the notation 

_ l r V®* — x v * + 


H v* ( V \ _ 1 f K <0* — x v * 4 CO Sin a ( B v 4 xC\ rnq 

« r . w ~si_l — <>“>■» 


(5 


from Eq. (5-7), it is easy to obtain the formulas 


0) F ~ + L 

v=i 


n gCXv—B^V ©» — x * » 




—<WTSs7ffx (5 

X [ (nr. jfg ~ *»* ± ” sin a , „ *v + M cos a ) ( D > + *C\sin 

X K « + v — » — M 1 - TvT"y C osO“)« to 

/*£- 2 <» + r* 5“ W - ^2 T v x ( 5 


v=l 


V = 1 


x t (x v 4 jo) cos a) f 1 — ) cos (ua) rfa 

' D V C / sin 

/ 3 A * hC—XB v H. 


‘ t ' 4 JXC 

dF, in / dh ~~ 2 Ll Yq jy (X v ) 


H-* ‘ 


v — 1 


(5 


a/ V<= / 5a) ^ C (to* -jr g y a* — x v ») -f j B u x v ,ff 11( . J. r 

= 2j ^ 


/ d(D 


j( 


<dV c 




oo — 2 


v=i - -V* — 2 nVc X (5 


X \ (gC 

— Tt 


x v + (ocosa V«*_x,*4 


2 4 <■> sin a ) / j i? v 4 xC^ C os 

® J \ 1 ~~DVC / sin ^ da 


We transform the algebraic factors in these equations. 
2 

A^C-B^ =0, with Eq . (5.21) taken into consideration, can be 
as 

(* + ~ x ^ 2 (i + i J ) + (H ix v ) 2 =0 ( 5 

According to Eq. (5.21) and (5.27), 

he xd v = h 4 ix v , g 4 V o ) 1 — * v » = (y^l— p ^~ x< . ( 5 


. 22 ) 

• 23) 

.24) 

.25) 

. 26 ) 

Equality 

presented 

27) 

28 ) 
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By applying Eq. (5.21) and (5.28)', 

gCy, - 2 ? v V<»'- x v » = (A + *x v ) |xg-(/31) v (5.29) 

can be obtained and, according to equations of the ( 5 . 12 ) and ( 5 . 28 ) 
type, 

c (“* + * V®*-V) 4 - s v x y = - (h + Xx v ) {h 4. (KZl) ¥ , yr+jj} i (5.30) 


Further, 
different , 


for P(K) = ic n +a 1 K n 1 +. 


. +a 


n- 


if all roots of k are 

v 


s 


x 


n-k 

v 


P’ (x v ) 


(1 (* = 1) 

\0 (A = 2, 3 n) 


The use of this theorem for n=4 with account taken of determina- /97 
tion of T v (5.18) simplifies the second sums of Eq . (5 . 23 )-( 5 . 26 ) ; 
only terms with factor (5-30) remain in Eq. (5-23), and only terms 
with factor (5.29) remain in Eq. ( 5 . 26 ). Therefore, with account 
taken of Eq. (5-28) with the notations 


N„ = BJ Vl+V, G,=,(h± ix v ) T\, 


(5.3D 


the formulas 


(o F 

•‘in 


= *2 


G v r Kh — x v n H 8 

o)/c [ x<r ~ {}r ~ 1)v y c J v (Xv) ± 


± 1 3in (ua) *L 

£) y Q J cos 2 ji 


/3ft ^ Cy //^ 

3/t, ldh-2 . J /c // (x «' ) 

W = 1 


/ 3? 

3^, / dg 


H, 


-“2 ^ 


v=l 


(Xy) 


j 3 /’^ / 3 Cl) 

‘dF^/d<o 



(Xy) “ 


1 




a v 4- X C 

z> yc 



are obtained and then, with account taken of the complex conjugate 
nature of the roots, 
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'£ £'-{*"■£ w}± 

* 

± 2g 5 cos ^ a) £& 2 lm w. (* v + * rrn?)> Oi * o) 

— 71 V— 1 

d^c / dh _ __ 2 _ J, r a* 

of^/dh yrr*^ Re \ (Xv) j ! 

dF v I'dg J, r % 

«W ** = “ 2 2 Im { c v^ (Xv)} j 

aF » c = h_ dF^ c /Bh_j_ BF^ C / dg i 

dF v.» / 30) o) dF^ U dh to BF^ I dg + 

2 71 2 ■ 

+ co/rrx* S If &“> 2 S 2 im (c v /v v (w v -f x v i + x*)} (n y= o) j 


dFo 


= —A dF o — £ gfo _ 


2 Im {CyVj 4- 


du a) 5A w dg co Y 1 + X 8 ^T x 

+ 0) y\ + X 2 S 2irD. 2 Im {C^v (IV. + '^H A. 2 )} 

— 71 V = 1 


(5.32) 

(5.33) 
(5.3*0 

(5.35) 

(5.36) 


We will therefore study roots k and functions of 

v 

v=l and v=2. 


k only for 
v 


Let a and b be real numbers. Then, 

V7+Ib = ^ 7 |(K V7* 4 6 2 4 ^ liKla’ + i’-a) (&> 0 ) 

Based on this, with the notation 

a£ = h { ^ - «* <* 4- **)]* 4 (1 4 X*) ± [<B» 4 A J - g 1 (1 4 X*)]) v ' (5.37 ) 

Y o>* 4 (h 4 i? /l 4 **)* = Ml 4 *Af a , A/iAf, = gh V 1 4 X 5 ^ 5 " ^ ^ 


Equations (5.13) and (5-16) can now be presented as 


*» = X (A 4 i* V 1 4 **) ± Sgn X /l 4 X 2 (A/, 4 *Af,) 

X J 

V 2 0)1 

j — V 1 + *■* + ig Vi + ^ 2 ) ± I A, | (A/ t + M/ 2 ) 

K x a * — <o a 


(5.39) 

(5-*»0) 


and Eq. (5-31) 5 with Eq. (5.21) and (5*19) taken into account, take 
the form 
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= — &g T sgn X (A/j + iM t ) 


G = 


(— l) v ‘ sgn X 


(— l) v sgn X (Af, + tA/i) 


v " „ \ — 0511 />. -r u fi V 

4 Ko)» + (A + ig V 1 + X 1 )* ~~ 4 V l®* + h*— T' (14- X»))> + 4AV d + X 3 ) 


(5- 4l) 
(5.42) 


By assuming e 


ia 


=C and applying the method of subtractions 


5 


JL [ cos ya da _ 1 f ? ft _ (e Vy* - 1 + y) 11 
2n J n COsa — V n ‘H 3 — 2j/S + l &Yyi~i 


(A = 0,1,2. . . 0!(5.43) 


can be found, where e=+l, in which the sign is selected from the con- 
dition | e*^=l+y | <l. 

Since 2 sin ya sin a =cos ( y-1 ) a-cos ( y+1 ) a, 

1 f sin yet sma^ = f° (^ = 0) (5.44) 

2n J cosa-y (-(eVy ! -l+# (y = 1, 2, . . .) 


Consequently, the first part of integral (5.22) will be 


where 


TT 
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1 /V — to 3 4 - cd sin a COs 
u cos a — x v sin 


(ya) rfa = 


ie 

0 


±‘1 


(y = 0) 

0* = 1, 2 ) 


<I> (xj = 


cV'xJ 


O) 


O ) 3 4- x 
1 = 0 


(5-45) 

(5.46) 


and e is selected from the condition |$(k )|<1. According to Eq . 
(5-39) and (5.40), 

(*‘ ) = (e V 1 + + X) [A 4- ig V\ +T\± e sgn X (A/, -f LV*) ] (5-47) 


Further, expression (5-6) can be presented as 


D i =a> i — 2u)L cos (a — a) 4- 4 - (A + Xr )* 

£ — y x * g 2 > cos a=xlL, sin a = — g/L 


(5.48) 
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Let 0=a -a. Then, 


D * = <■>* — cos Q + IS + (h + X x )t 

cos (fia) rfa = ? - x„« + m sin (0 + a) COs 

( cos a xj /> sin j [co cos (0 + a) — v ) /) 


sin 


(5.49) 


~f o) I/O] 


It can be proved in succession that 

1 o) cos o cos ft — x Hk co sin a sin ft 

co cos(ft + 0) — x (co cos 0 )* — 2 x cos o (co cos ft) x* — co* sin* o j 

V co * x a + co sin (ft 4 - a) _ YjSj^ x 1 cos 0 — x sin 0 + co s in ft 
co cos (ft + o) — x (a cog ^ x cos <3 — a)* — x 1 sin a ~ 

On this basis, the second part of expression (5.22) will be 


jAo*— x v » + co sin a c03 
2 n (co cos a — x v )D sin ~ 


«=f; 


Tk 

co r 

;<w«S 


sin ^ft sin ft rfft 


cos }/ co 2 — x v * cos a — x v sin o 

(aa) 

sin co 


J (co cos ft — x v cos a — jAo a — x v a sin a) 

Tl 

^ COS flft (fft 


+ 


J (fa) COS 0 — x v cos a — }/~co 2 — x,, 3 sin a) D 
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As a result, after substitution of 8 = tt— 2«J> , 


H 0.. v = 0, J* (x v ) = ± “ (X v ) - (- i)“ El ^ + Xa) : 

1 « ^(f> + t) 3 + (A + Xx)» 

f cos — X v 2 cos a — x v sin 0 oj n A ^ 

X 1 Ow)- 77 n „ (*. Pv) ± co3 (pa) (Ac, p v ) } 


is obtained, where 


** =r 


4 (oL 


2 co 


(o>+L)*+(A + **)* ’ ^=- u4 . XwCOSO+ 


rt/2 

(Ac, p) = ^ - 


cos 2pq> rfy 


5 (1 + A>sin» <p) Kl — A 1 sin*9 ’ ^ ^ 


(5.5D 


sin a 

Vi(*,P)- (fc,p) (5.52) 


Function n^(k, p) is a particular case of the generalized el- 
liptical integral introduced and studied by the author. In the Ap- 
pendix, recurrent formula 

'A {ry , (Ac, P ) + ry x (Ac, p)} = (i + 2/p) (*, p) - (2 / p) (*) (5.53) 


33 


was obtained and convenient formulas for calculation of ^(k, P) are 
given for any (real and complex) value of parameter p. 


Further, based on Eq. (5-^8) and with the use of the substitu- 


tion a-a=TT-2<}), we obtain with module k (5-51) 

n 

da 


cos 

sin 


<|W> 


2 


2nD n V (<d + L) 1 + (A + X*) 1 sin 


(5.5*0 


In the case of g=0, when M 1 =v / to 2 +h^ and M 2 = 0 

O, = sgn X (VT+~xHM), s = — sgn X 

y b) a + h* + h (y t + >.i — I X|), e = — sgn X, ( h 1 < oW) 
0) 

_ sen X V^+h' — h (/ITT 1 + 1X1), e = sgn X, (A J > oi s X>) 


f.;:V 

d>i 

W.',’ - . 
ri*l 


‘ 


v i k . 

— sgn 

6,= . 


y 

— sgn 


(5.55) 

(5-56) 


(x u ) = ieO/ + i 


)/x v ‘ _ o)* (_ sgn x)* +l Py (N v + x V 1+ >.’) 


<0 


n ^(to + L)» + (h + Xi)» 


n„ (x, p v ) 


// , w _ _ <*, = o (5.57) 

^ Jt / (CO + L ) 1 + (fc + X*)* * 




4(oL 


Pi 


= =F 


2co sgn X sgn r 


(CO Hr £)* + (* + M* *> a ^ /co* + h 1 V 1 +X a ± 1 1 \h \+ cosgnXsgnxl 


and with the notations 


Si o r / 1 -f X 1 £F sgn X 'Kco 1 /t» 

£a l^co 1 4- A* V 1-j- X* i [ | XA | + cosgnXsgn x] (5*58) 

z z \ = a vt+t> ± rx i 


the following final formulas are obtained from Eq. (5-32)— (5-36) : 


HF e = + Of _ (~ sgn g) tttl sgn X (X. Pl ) - (X, p,) 

2 Vl + X* Vl + X* “ n v (0) + L) 1 + (A + Xr)i 

dF ye = _ 0)^ + sgn (fti — QiX 2 ) O/ ^ 

dh 2 V 1 + X* V«>* + A* + 

4 . (— sgn x)^ z^n^x, p t ) + (x, p,) 

Vl + X* V toi + A J n V (CO + L)1 + (X + Xx) 1 
=-. r i , (- sgn x)^ 1 S t ny (X, p,) + g a n„ A (X, p a) 

dg .2 Vw 1 -fA* Ko> l +A J n V'(co + L)> + (A -f Xx)» 


(5.59) 
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3*4 


dF 0 

d& 


d Jj ± C = _ h dF„ . 2 (- sgn x)^ x/to ^ (->/,) (Jt) 

ato a dh n ^<o4-L)*+ (/»4-Xx) 3 11 

_ A + JL g/ <i) r.(-v.) (W _ — . 1 ■— 

to dh n f(w+L)>+(H Xx) 3 to v 1 + x* 


The remaining functions are absent. Since parameters P v (5-57) 
satisfy the condition l+p v >0, the corresponding integrals of the third 
kind of n(k, p ) can be expressed through integrals of F(k, a) and 
E(k, a) by Eq. (6.1)— (6.4) specified below. 


If x=0 , and, in Eq . (5.22), 

b v + xC N y + x VTTw 

~©TF 5 


1 — (— l) v sgn X 


There are obtained as a result the simple formulas 


2 V l + = (1 -f sgn X) O/ + (1 — sgn X)^ 


2 K 1 + X 3 "Kto* + = — 

dh 

r _ dF 

sgn k 2 / co 2 4 - h 2 — v±_ = 
dg 


dF 


\LC 


h dF ] 


\LC 


did 


(1 + sgn X) tD/ — sgn (h 2 — to s X 2 ) (1 — sgn XJ02 11 

— (1 + sgn XJO^ + (1 — sgn X)®, 11 

a^ = _ h dF 0 l 

ato to dh u Y 1 4. X 3 


( 5 . 60 ) 


The remaining functions are absent. 

Equations (5-60) are simplified in the important case of x=0=h 
with A>0, 

Ot 11 dF^ , 

“aF = TT’ o,= Vi + i?-x 

d J\* _ J — 1 /to v 1 4. jj (n = 0) ( 5 . 61 ) 

a<o I 0 (i» = i. 2, . . .) 

In considering the case of x=0=h, it must be noted that it is 
evident from Eq. (5-2) that, for h=0, functions , 3F^/3g, 8F^/3to 
and A3F^/3h do not depend on the sign of A. Functions only of A>0 
will therefore be considered subsequently. 
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According to Eq. (5-48), for x = 0, L=g, cos a=0, sin a=-l; there- 


fore , 


= (- 0“. 


cos COS / Jt \ 

sin (na) = (- 1 ) 1 ‘ sin (pt) 




k* = 


Hos ( x v) — 0» u (x v ) = i \ G** 1 (x v ) — 

X1 \kS A 

{r (•* j ) it "► <*• '•> ± % (i*f ) v <*■ *> } 

4c aL 


(to £)* 




2co 


2co 


(0 — ^ G) a — x v 1 w — l* V X y 2 — CO* j 


( 5 . 62 ) 

(5.63) 


It is useful here to use the transformation established by the 
author (see Appendix) 

l + k* p» ( ( k, 2 \ fc, 1 

n,(i lPv ) = — i--) n 0 (A:,.,,) - — /•(*♦)} 


*• = 


k* 


1 — 1^1 — Jfc a = _ 

i + Vi — *> (i + 


i . p* = 


1^+^ T = _( 1+ A) (5 .6iO 


According to Eq . (5-63), 

j/co’-V 


T = 




fco /L< 1 
* = \l/co<1 


(5.65) 


and, for g=L and h=0, with Eq. (5.27) taken into account, the result 
is 

_ ^ 2£. y p) 1 — x v a + L 3 co s / L’ 


/>* = 


x v 3 L* / <o 3 

co> -f 2L y <o*— x v > 4- Z, 3 


1 + i 3 

1 

i + X* 


(co<L) 

(u>L) 


( 5 . 66 ) 


Parameter p* turns out to be real and the same for all k . 
Moreover, for x=0=h, function F and its derivatives have two branches 

v 

each 

(a) £ J (1 4- X 3 ) a 1 > 0, (b^ 3 (1 4- X 3 ) -a) a < 0 (5.67) 

where quantities < v and p v are complex. in case (b). An extremely 
simple case of change to integrals of the third kind with a real pa- 
rameter thus occurs here (see Appendix). 


36 


In case (a) 


2g) 


r;Afi = 0, M t = V & (i -p X*) — o>» * ^ “ to 4-Z, (1-f X*) ± X A/, 

; «©(*)=£ (e /r+T* + X ) (L V r+T» ± eA/ s ) 


( 5 . 68 ) 


Since lA+A 2 = v^^+co 2 , it follows from condition |$(k )|<1 that 


® (xj = — £ L V^i -f x»— M± (y i + x j — x) 
0) 


8 = — 1 


«*(**)= j 


. L V 1 + X 1 + Afi (y < _p x* — xh e = — 1 (<oX>Afj) i ( 5 . 69 ) 
0) 

. L V 1 + X 1 — Mt (yj .4. X» 4. x). e = 1 (o> X < Af*) 


2 2 2 2 2 2 

However, Mp -to A =(L -to ) ( 1+ A )• Consequently, 


O (x.) = £ej ® (x v ) I, 


e v = — sgn [w — (— 1) V L] 


(5.70) 


and, for L>to , the first branch of $(< 2 ) is missing. Parameters 
(5-68) are real, and transformation formula (5-64) contains only real 
quantities. We designate 

Tx _ 2_ y\ + X» (M , ± XL) 1 

(5*71) 


Ti n a>+L o> 4 - L (1 4- X J ) ± X Mt 
5i 2 0 MjdbXX ai =A/,±XL 


6 , n g> + L w + L (i 4 X«) ± X Af, 




Prom Eq . (5-63), it is easy to obtain 

?M C ) ■ - (< TI •( : ) r + « (rf) 5 (*• S)T «*> 0“r) t "/ (*• ?.) 

o-(« T •( : ) r + -(»-f ) ; (*• s ) ± ‘-k-) £ ->* (*. s 


The final formulas will be: for y = l, 3, 5, • • . , 


a/ 1 ,. 


^ c = - XL_£L+(-l) a 


^ f aiTin^ (fci.pi) + a»Ta 11^ (fc.Pi) 2 £, 


I 2 1^1 4 X s VTb’ (1 + X S )- u 1 


— — £ (-•/•) /A 

j; to-pL^V ' A> | 


* dr + { i} 2 vr+vVL'a+\>)-o* 


dF. 


Jl+I 


VAC 


dh 

dF, 


_ ( i) 2 — 1 Oil* 1, + 1 O2 P 


2/1 + V Vl*(\ + W) — o* 


(5.72) 


VA + 1 


^ = n ~ Ti^ (A, Pi) — T»ry (Ac, pa) 
a/l 2 /m* / L* (1 -Jr M) — o>* 


_ r t (*■ fO + (x, a/y, exiot |^- e> |®x|' 

d * 2 /L* (1 + X*) — 0 )* 2 V L> (1 4 x j ) — <o> 


11+1 
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and for y=0, 2, , . . 


»F ur J — «i|q>ir^«»iq>,r 

dh ’ 2/i ^X,>/Z.*(l+X»)-w> 

tt F = i f. d Jjil U - f aiTlI V (k ' Pl) + a ir * n u (*. Pi) 2 L ) 

5/ * (2 Kl 4- X^Z^i+X 2 ) — <o 2 n(Hi) V' 7 ’ Wj 


ai | <D t I 1 * 4- ai | P 2 | n 




dF 

V-c 

dh 

OF 

t*a 


= < ~ TlI V (k, P i) + Tin a (k.pi) 

2 V 1 + x» yi^l + X 2 ) - o>» 


ax 


= (- 1) 2 


ap 


I p# - 1 0,1 * 

2 y 1 + x* yz* (i + - to*' 


af = ( - 1 > 2 


| <Pi | **• — e, | < p, \v- 
2^(1+ X 1 ) — to* 


fAi = f _ (At, + fcn/p, Pl > 

^ Z V 1*11+ X 2 ) — 0)> 


>¥=0) / (5.73) 


Moreover , 


’a/^ / a<o 

aFjj, / a&> 


a/^ / a? 

<0 a/^, 1 d g & = 1. 2, 3, . . .), 


dfj 

cHd 


L dF„ 1_ 

“ dg co Y 1~ X* 


(5-74) 


In case (b). 


Mi = V w» — L* (1 + X 2 ), Aft = 0 

Xl =(u ■£> dt f a^i) yr+T> = - yr+T 2 ^ 


*3 


A r 2 


■ *rnPS no{«. ",r w}: 

v=l V / 

±- Sin ( 11 JL) — — j — £ (-7.) m 

X n cos y 2 j «> + ^ !*• W 

v, / aa yr+l? ^ 2 M X lm (// (O f 


dF^/dg 

aF V.> 


a^* c _/ao) 


^ ds _ vi ( b v n J Uy-c . . 1 
/aff 2 j 2a/ ! *«{* K) 

v*i ^ y ^ 


a^^aa,; 


<» a/’n, / a* • 


ap. 


L dF„ 


“ °S coV l 4. X 2 

“ <1> (x!)= (yi + **-*•) (± A/ t - iL V 1 + X») 

e= — 1, 1 0> (x v ) I = y 1 + X 2 — X 

2o) 


Pi 

Pi <0 + L (1 + X 2 ) ^ iX A/ t 

2 “ - fo-f Z.(l + X 2 )l =F & Me 


K“ + Ly 


1 + x* 


= p (1) i iP <2) 


(5.75) 


(5.76) 


(5.77) 


(5.78) 
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It is difficult however to separate the real and imaginary parts 
of functions H^ c and H^ s in general form (for any number y) in Eq. 
(5*63) since, for a specific number y, it is initially necessary to 
express n (k, p v ) through n Q (k, p^) by means of Eq. (5-53) and to 
use transformation formula (5*64) with complex p . 

6. Appendix 


In the preceding text, complete elliptical integrals of the 
third kind of n(k, p) are found, parameter p of which is complex or 
changes from -1 to °°. If p>-l, II(k, p) can be expressed, based on 
the Legendre equation [4, pp. 402-403], by elliptical integrals of 
the first and second kinds of F(k, 3) and E(k,B). With the notations 


&4 ; Qi (*, P)=y-/(«B (*', p) - E ( k ) F{k',P) + F (k) F (k\ P) ( 6 1 ) 
Qi (*. P) = F (*) E ( k , P) - E (k) F (k, P) (*' = V 1 - Jt 1 ) 


JPiV*- • ■ 

Si : 


«JT/, v I/ ~ Ql (*. P) . f 
|D(*. P ) = |/ , +p? == + T 


.... r , , F(k) 1 I 

+ T+i- ,mp “7T+7 (p> > <6 - 2) 




{i/~=T Mh P) 

! V t+~p i n 

Qi (fc. P) 


+7/i5+7 sin|J = Xzi£ «><-,<*■> (6 _ 3) 


l 


Vt 


! 1/ — p 


+ pV— (p + k*) 


F(k), sinp=]/i±il (*’<-/><!) 


According to Eq. 
will be 


(3-45), in the intermediate case (p=-k ), 


n ( k , 



rftp 

(1 — k 2 sin J 


E (ft) 

1 — ft 1 


(6.4) 


it 


It is easy to see that (^(O, g) = 0, Q ? (l, B) =co , 


Qi (*. P) 


*/,«- P (*=1) 

»/,«(! — sin P) ( ft = 0 ) 


0 < Q, (ft, p) < n / 2 (6.5) 


The case of a complex p can be reduced to the case of a real one 
by means of the transformation found by the author. The results will 
be presented briefly here. 
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( 6 . 6 ) 


For transformations, in place of integral 



<*£ 

(1 + p sin* <p) V 1 — k * sin* <p 


(** < 1) 


it is more convenient to consider integral 


J 



c[9 

(e — sin*<p) Vl — k 2 sin*<p 


— n (Ar, — e a) 

S 


(6.7) 


The derivative with respect to <{> of function 

( 6 . 8 ) 

(u<0) 

can be presented in this form 

d .k _ (. Q (J) ) l 

d <? X k*R (t) ) Vi — k 2 t ' 1 ~ Sln J<P 

where Q(t) is a polynomial of the second degree, 

R(t) = t* + 2k '” 1 t » + t - IL (6-9) 

k* k 2 k 2 

With the corresponding numbers A, B and E, 

R (0 = (<* + At + B) (t — E) = t* + (.A — E) i* -f (B — AE) t — BE (6.10) 


1 


T. ;<r) = 


Vu 


“ arctg 


sin 9 cosqp V u 


(u>0) 


[2 V-u 


In 


(p + sin a <p) V 1 — /c*sin a 9 
V — u sin 9 cos 9 + (p + sin* 9) V 1 — k 2 sin 2 9 
V — u sin 9 cos 9 — (p + sin*9) V 1 — k 2 sin 2 9 


The roots of polynomial R(t) will be 


e i A 

e, =-~T±iVs, 




e, = E 


( 6 . 11 ) 


and it is easy to establish such a connection between three inte- 
grals of the Eq. (6.7) type, 

k 2 L (a) = k 2 F ( k< a) + IF/ j (e v , a ), w = (6.12) 

v = l R’ (e v ) 


If numbers A and B are fixed, number E and parameters u and v 
of function (6.8) are subject to determination. Comparison of Eq. 
(6.9) and (6.10) gives the three necessary equations and exclusion 


\ 


of u and E leads to a quadratic equation for v. As a result. 


»i 

Vi 


(1 - A») B ± VK 
1 + A + k'B 


k*E = 

D 


u = k* (A — E) — 2k'v 
K = £ (1 + -f D) (1 + A l l + k*B) 


( 6 . 13 ) 

(6.14) 


Let roots and e 2 (6.11) always be complex; in this case, 
S>0, and it can be shown that K>0 and the numbers of Eq. (6.13) are 
real, E-^E^l, 

Wi= Q(E)/R'(E), 

R'(E) — B + AE -f- £*. Q (E) = (2 — E) k*E* — E 4- t> (1 — 2E + A*£*) 


m 


M ± iN, 


M = 


<r_ 4- *p 
^-fp 1 ’ 


w _ xr — ? p 


(6.16) 


<7 = V a .4 [1 - AIB 4- (2 4- 4) 4- * (2 + A14)] — (2 + yl) A»fi -f p (1 — A*£) 

x = — [1 — k?B + (2 4- A) k*A + v (2 -f *t4)] VS (6.17) 

r = — 2S, p = — (A + 2E) r* 4. p» = (£) 


Further, with S>0, 

J (e lt a) = / (1) (a) + U (2) (a), J (e„ a) = J M (a) — £7 <2) (a) ( 6 . 18 ) 

and the real parts of Eq . (6.12) form the equation 

MJ m (a) — NjW (a) = »/, {k 2 L (a) — A *F (A, a) — \V 3 J (E, a)} = U (a) (6.19) 


By taking two roots v 1 and v 2 (6.13) and numbers M, N and U 
corresponding to them, two equations of the (6.19) type can be 
written. The solution of such a system of equations will be 

J M = (- U t N t + U t Ni) / A, J (2) = (U 2 M X - U X M,) I A, A = - + M I N 1 (6.20) 

Since the integral of J(e, a) of the (6.7) type with complex 

2 

e, analytical function e or A and B (for k <1 even with a=Tr/2), 
quantities and exist in the case of A=0 (if it is possible) 

p 

In any case, with A=0 (as well as with 1+A+k B=0), one of the subse- 
quently indicated transformations of J(e, a), of which the first 
changes only e, can be used first. 

Thus, for e=a+ib, by finding the numbers 


41 


/105 


A = — 2a, 


B = a 1 4. b* 


( 6 . 21 ) 


J(e, a) can be presented in the form of a linear combination of in- 
tegrals J'(E^) and J(E 2 ) with real parameters E 1 >E 2 >1. In the case of 
complete integrals (a=ir/2), functions L-^(a) and L^(a) are absent, and 
integrals J(E) can be calculated simply by the above specified equa- 
tions (6.1)- (6. 4). 


I 
i 
! 

' ( 6 . 22 ) 

(X>0) 

it can be determined that 

d6 1 — 2 tA- h*t l . . v 

= r ~ (t =r Sin 3 9) 

d<P 1 — (*» -f X) t 4- Xi» Vi — k*t 


By differentiating over <f> the function 

arc tg (X<0) 

Y— X Y 1— Ac 1 sin 2 9 

0 (X, q>)= 

, 1 K ^ sin 9 cos <p + V 1 — ft 3 sin 3 9 

k 2 Y^ Y ^ sin 9 cos 9 — 1^1 — k 7 sin 2 9 


and the equation 

4 0 a ) = <*• “) + (1 - W J («. <*) + (!- **& ) / (6.a) 


(6.23) 


can be reached, where a and b are connected by the relationships 

1 — a 1 — 6 


b = 


1 — frV 


1 — k*b 


(6.24) 


From this, with notations p=-l/a, p#--l/b, it follows that 


pp*0 (pp», a) = k*P ( k , a) - (p 4- * 2 ) n (Ar, p, a) — (p* + fc ! ) n (Ar, p*. a) 

p + ft* p* + P 

P* ~ p + 1 ’ ^ — p* + 1 


(6.25) 

( 6 . 26 ) 


We call this transformation symmetrical linear-fractional. For 
k 2 <-p<l, there will be 0<p x <co, and the transformation proves to be 
extremely useful with p^j-1. With cx=tt/2, function 9 (6.22) is absent. 

We point out still another important transformation. We in- 
troduce variable 4 ^, with the use of the Landen substitution (pro- 
posed by him for transformation of elliptical integrals of the first 
and second kinds) 

42 
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tan;q), = - 9in 2 9 

*1 + cos 29, ’ 

dy 


hi 


= 1 — yi — a* 


a* 


1 + v 1 — ** (i + yi — **)* ’ 
L±A <*p t 1±A 


(1 + *,)> 

1 


/1-A»sin*q> 2 YY^~k£ sin* 9, ’ = FTT?^ 


By converting the first equality of (6.27) to a quadratic equa- 
tion for cos 2<f> , we obtain 

cos 2<p = — A x sin* <p t + cos (jij rr=r Vsin , 9i = l _ 23inS<p (6.28) 

The root was selected here which corresponds to 0<4><7r/2. Ac- 
tually, ^ $ 3 4>j = 0 with <J)=0, <P^ = tt with <j>=tt/ 2 . It can now be 

shown that 


a 

S 


^9 l+Ai IT cos q> t - dpi 

(e — sin* q>) V 1 — A* sin* 9 ~ (1 + Aj)* — 4A ie U e, — sin* «p 1 


*» a t 

+ (l-2e + e l A,)\ A* — k \ - rf 9i 1 (6.29) 

J (e x — sm*^,) Y\ — Aj* sin*^ 1 J Vl — A x * sin*^ J 

4e (1 — e) e (1 — e) .. s ; n 2a 

e ‘ - (1 + A,)> - 4A t e - 1 _ A*e <* + V 1 ~ **)*• tg °i = ^T + cos 2a (6.30) 


Let a, b be real numbers, e=a+ib; parameter e 1 will be real only 
under the condition 

(1 - 2a) (1 - A*a) -f A* [a (1 - a) 4 J*] = 0 


By changing to integral (6.6), the following relationships can 
be obtained: 

n (A. p, a) = L+A {(1 + A - A) n (A lt Pl , a l)+ Af (A lt aj + (6.31) 


+ 


c os 9 ,^ 1 * I-/*-** Pi = 1 - 2A lT + A,* /, , 2 \ 

J 1 4- Pi sin a 9 X J 1 + 1^1 — h 2 T 1 — 1 ’ T \ 1 + "F j 


0 (k ' P ' a) 2 (1 + V\ - A*)’ {/, + ft* F ( * 1 ’ “ l) + 7p4^Hp + M*) n (/:it Gl) 


+ 


+ 


P + 


— L ( p a.)V p. = — g A ki 

1 J pil(l + V"l - A*) 2 


( 6 . 32 ) 


(o, A) = f cos _ 1 i 14- t Yq sin ft 

)l4?sin 2 ft ,2 Yq l-£/7sin0 


il3 


(6.33) 


Since k-^k, by using transformation (6.32) n times, an ex- 
tremely small module k n can be reached, and it can be assumed with 
great accuracy that 


n tv „ d<Pn 1 , 1 ± { V 1 +Pntgq n 

n. Pn> n) J 1 + Pn si n’<p n f2 yT+J^ ln 1 - i yV+T n tg a n 




p 

Modules k decrease extremely rapidly (for example, k 0 <10 

2 n 2 -8 2 ^ 
for k =0.6 and k^ <10 for k =0.99) , and the successive approxima- 
tions method proves to be quite effective. 


For complete integrals (a=ir/2), we have 


a n = 2 a„_ lt L(p„, ji) = 0. 2 F (k) = n (1 + *,) (1 + k t ) . . . 

P 3 + 2 P + * 3 




(p + i) (p + fc 3 ) 


n (k lt p t ) 


} 


(6.35) 


It can be shown that, for real numbers g and h. 


M = ( ‘/ a « (l + gY' 1 ' (i + g > 0) 

1 -f g sin 3 0 \ o (1 + g<0) 

dQ Nj — t sgn h .V, ji 

1 + (g + ih) sin 3 0 VW+ gY + * 3 2 


( h =h 0) 


_ f/(l + g) 3 + /< 3 ±(1 If) l 7 ’ 
N 2 \ 2 J 


In generalizing elliptical integrals of all kinds, the follow- 
ing function should be introduced, 

n m <n) (t, p, a) = [ £ ° s2 "y (i _ k * sin i (6.36) 

.) i “h P sin 2 9 
o 

Since 

cos (2m-2)(p + cos (2m + 2) <p = 2 cos 2m<p (1 — 2 sin 2 cp) (6.37) 

V, (n^”’ (k,p,a) -f nj"’ (At, p. a)} = (1 + 2 / p) n m (n) (A •, p, a) — (2 / p) /T^ n) (k, a) 

Em^ n ^(k, a) is a function introduced above (in Section 3)> 
generalizes integrals of the first and second kinds. 


then , 

where 

which 


/107 


REFERENCES 


1. Maykapar, G.I., "Applications of vortex theory of the screw," 

Tr. TsAGI 613, (1947). 

2. Proskuryakov, A.P., "Theory of the rotor at zero angle of at- 

tack," PMM 20/4, (1956). 

3- Bratukhin, I.P., Avtozhiry (Teoriya i raschet) [Autogyros (Theory 
and Calculation)], Gosmasmetizdat Press, 1934. 

4. Gradshteyn, I.S. and I.M. Ryzhik, Tablitsy integralov, symm, 

ryadov i proizvedeniy [Tables of Integrals, Sums, Series and 
Products], Fizmatgiz Press, 1962. 

5- Zhuravskiy, A.M., Spravochnik po ellinticheskim funktsiyam 

[Handbook of Elliptical Functions], Moscow-Leningrad , USSR 
Academy of Sciences Press, 1941. 

6. Samoylova-Yakhontova, N.S., Tablitsy ellipticheskikh integralov 
[Tables of Elliptical Integrals], ONTI , 1935. 

7- Belyakov, V.M., R.I. Kravtsova and M.G. Rappoport, Tablitsy 

ellipticheskikh integralov [Tables of Elliptical Integrals ] , 
USSR Academy of Sciences Press, Vol. 1, 1962. 




